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We study the phenomenon of intermittency in an inhomogeneous lattice 
of coupled maps where the inhomogeneity appears in the form of different 
values of the map parameter at adjacent sites. This system exhibits spatio- 
temporal intermittency over substantial regions of parameter space. We also 
observed the unusual feature of purely spatial intermittency accompanied by 
temporal periodicity in some regions of the parameter space. The intermit- 
tency appears as a result of bifurcations of co-dimension two in such systems. 
We identify the types of bifurcations which lead to spatio-temporal or purely 
spatial intermittency and identify the mechanism of intermittency for one 
such case, that arising from the bifurcation of the synchronised fixed point 
solution, by examining the spatial return map. The intermittency near the 
bifurcation points is associated with power-law distributions for the laminar 
lengths. The scaling laws for the laminar length distributions are obtained. 
Three distinct types of intermittency characterised by power-laws with three 
exponents which fall in three distinct ranges can be seen. The structure associ- 
ated with the tangent-period-doubling bifurcation undergoes further bifurca- 
tion to a spatio-temporally intermittent regime with an associated power-law 
exponent. Three of the exponents seen in this model show very good agree- 
ment with those observed in fluid experiments involving quasi-one dimensional 
geometries. 
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The existence of spatio-temporal as well as temporal intermittency is an interesting phe- 
nomenon which is found in a wide range of dynamical systems of physical interest such as 
oscillators 0, chemical reactions 0, pattern formation of various kinds 0, turbulence ||, 
fluid flows lH and signals of all kinds. The phenonomenon of temporal intermittency has 
been studied extensively and is relatively well understood §. Pomeau and Manneville @ 
first established the route to chaos via temporal intermittency and studied its behaviour. 
The nature of intermittency in spatially extended dynamical systems, however has not been 
understood very well. The presence of spatial as well as temporal intermittency has im- 
plications for understanding the physics of pattern formation and for understanding the 
ubiquitous presence of structures in chaotic systems. 

Coupled map lattices are particularly simple paradigms for systems with extended spatial 
dimension which show a wide range of interesting behaviour ||. Spatio-temporal intermit- 
tency, which has been defined as a 'fluctuating mixture of regular and turbulent domains 
U', has been observed in such systems M || | lQfl . However, the mechanism by which such 
intermittency occurs has not been very clear. We study the phenomenon of spatio-temporal 
intermittency in an inhomogeneous lattice of coupled maps where the inhomogeneity appears 
in the form of different values of map parameter at different sites. Such lattices have been 
considered in the case of pinning studies [11]] and in the context of control of spatio-temporal 
chaos ||10|| . The coupling between maps is diffusive and nearest neighbour. We consider the 
simplest case, that is, a situation where adjacent lattice sites have different values of the 
parameter and alternate sites have the same value of the parameter. 

Our model is defined by the evolution equations 
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where f(x) = — x) is the logistic map and /i, /i'e [0, 4], x\ is the value of the variable 
x at the lattice site i at time t, and < x < 1. We set // = /i — 7 and use periodic boundary 
conditions where 2N, the number of lattice sites is even. The synchronised fixed points of 
the system are given by x = and x — A '~ 7£ " 
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Expanding the evolution equations about the synchronised fixed point, the linear stability 
matrix has the form 
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This matrix is of the block circulant form. Thus it can be put in a block diagonal form 
by a similarity transformation [|T!J, 13] of blocks M{6) where M{6) are 2x2 matrices of 
the form 



M(0) 



(l-e)f'(x) Ul + e^)f',(x) 



;i + e- 4e )/:(x) 



MM*) 
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where, 8 



h...N. 



2n(l-l) j 
N ' 

The eigenvalues of this matrix as a function of 6 are given by 



(1 " + #0 ± ((1 - + f'A 2 ~ 2(1 - e )e(/;/;0 2 (l + «w(6>))) 
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(4) 
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It is clear from Eq. |] that the eigenvalues of the matrices M{9) are bounded between 
the largest eigen- value of the matrix M(0) and the smallest eigen- value of the matrix M{n). 
For the x = fixed point, the eigenvalue equations are 

A 2 - A(l - e)(2/i - 7) + /i(l - 2e)0 - 7) = forM(0) 

A 2 - A(l - e)(2fi - 7) + Ml - e) 2 (/i - 7) = forM(n) (5) 

For the fixed point, x = ^~2^,T 1 the eigenvalue equations are 

A 2 - A(l - e)(2fi - 7)^ + Ml " 2e)(/i - 7 )£ 2 = /orM(O) 

A 2 _ A(l - e)(2fi - 7)^ + Ml - e) 2 e 2 (/i - 7) = forM(n) (6) 

where £ = 2£^2. 

We fix the value of [i at 4.0. The system now has two parameters e and 7. It is clear 
from Eq. H that the eig en- values of the linear stability matrix are functions of e and 7 . 
It is well-known [14| that bifurcations occur and the nature of the stable solution changes 



at the values of the parameters where the modulus of the eigen-values of the linear stability 
matrix crosses the unit circle. In the co-dimension one case, the tangent bifurcation can 
be seen where the largest real eigen- value crosses +1, the period-doubling bifurcation can 
be seen where the smallest real eigen- value crosses —1, and a Hopf bifurcation can be seen 
where the complex eigen-values cross the unit circle. For the present system, the conditions 
for the eigenvalues to cross the unit circle define a set of curves in the two-parameter e — 7 
space as listed in Table I. These curves intersect in several places where the equations 
are simultaneously satisfied. At these intersections two eigen-directions become unstable, 
resulting in bifurcations of co-dimension two. A rich variety of spatio-temporal behaviour 
can be seen in the neighbourhood of such points. We concentrate on regions which exhibit 
spatio-temporal intermittency where a fluctuating mixture of regular and irregular domains 
can be seen. Such solutions are possible over a large region of parameter space. Table 
I lists the possible bifurcations associated with intermittent behaviour and the different 
types of intermittent behaviour observed near such points. It is interesting to note that 
in addition to the phenomenon of spatio-temporal intermittency, the inhomogeneous lattice 
has regions in parameter space where pure spatial intermittency accompanied by temporally 
periodic behaviour can be seen i.e. the temporal behaviour is periodic, but spatially laminar 
and turbulent regions co-exist at any point in time (See Fig. 1). The temporally periodic 
spatial structure seen shows the unusual occurence of long range spatial correlations. This 
phenomenon can be seen in the vicinity of the tangent-period doubling bifurcation whereas 
all the other bifurcations result in spatio-temporal intermittency. The bifurcations in case 
of the Double Hopf and period-doubling are inverse bifurcations as can be confirmed by a 
normal form analysis. The details of this calculation can be found elsewhere [ |T9| . 

Table I shows that spatio-temporal as well as purely spatial intermittency can arise in 
the neighbourhood of co-dimension 2 bifurcations. We examine one such case, that of the 
tangent-period doubling bifurcation where the interesting case of purely spatial intermittency 
can be seen. Let Xo(7> e , A) be the characteristic polynomial for M(0) and Xtt(i, A) be the 
characteristic polynomial for M(n). Then the conditions to be simultaneously satisfied for 
the Tangent-Period doubling bifurcation are 

1 - (1 - e)(2// - 7 )£ + Ml - 2e)0u - 7 )£ 2 = 

1 + (1 - e)(2fi - 7 )£ + Ml " e ) V - 7)£ 2 = (7) 

where £ = 7<S ^ M ~ 2 as before. These conditions are simultaneously satisfied when 7 = 
1.18, e = 0.63. ^ ^ 

We show the spatial intermittency present in the vicinity of this point in Fig 1(a). The 
size of lattice chosen was 1000 and it was iterated for 10, 000 iterations after discarding 



3 



20, 000 transients. The ordinate is the value of x at the zth site and abscissa is the site 
label. The lattice is plotted at a fixed time. In Fig. 1(b) we plot the time evolution of the 
lattice which shows stable period two oscillations. Thus the intermittency is purely spatial 
in nature and temporally we have stable periodic behaviour. The temporally periodic spatial 
structure seen shows the unusual occurrence of long range spatial correl ations. This can be 
seen in the distribution of laminar lengths which we will discuss shortly. 

Another tool that can give some insight into the intermittency for this Tangent-Period- 
doubling case is the spatial return map plotted in Figure 3. Before the bifurcation, the return 
map shows only a single point corresponding to the synchronised fixed point solution. After 
the bifurcation a loop is seen. This is because intermediate values between the two period 
doubled points are also possible now. This loop acts as a relaminarisation mechanism and 
gives rise to intermittency. The region in which the intermittency occurs in inhomogeneous 
lattices is considerably bigger than that seen for homogeneous lattices. This is because the 
lattice is no longer symmetric under the "flip". Thus the periodic solutions which were 
stable for homogeneous lattices bifurcate and give rise to intermittency. 

Each type of bifurcation gives rise to intermittency of a distinct type as typified by the dis- 
tribution of laminar lengths. The distribution of laminar lengths shows power law behaviour, 
and the tangent, the period-doubling and the Hopf bifurcation are associated with their own 
distinct power laws. The distributions for the tangent-Hopf and tangent-period-doubling bi- 
furcations scale with the power-laws associated with the Hopf and period-doubling bifurca- 
tions respectively. These power-laws are seen in the vicinity of the co-dimension two points 
listed in Table I and also in the neighbourhood of the lines which correspond to the tangent 
bifurcations. The question of whether spatio-temporal intermittency persists beyond these 
regions and whether the distribution of laminar lengths in the spatio-temporally intermittent 
regime crosses over from power-law behaviour to exponential behaviour is presently under 
investigation. 

To distinguish between various kinds of intermittency, we calculate the distribution of 
laminar lengths. The length of the laminar bursts, i.e. the number of consecutive sites which 
follow periodic behaviour before being interrupted by chaotic bursts is calculated. The distri- 
bution for this length shows a power law behaviour with P(l) fh r. This shows the presence 
of long-range spatial correlations in the lattice. The power law exponent characterises the 
type of intermittency. The distribution of laminar lengths for three kinds of intermittency 
is shown in Fig. 2. The exponent Ci is associated with the tangent bifurcation,the exponent 
(2 with the Tangent-period-doubling bifurcation and period-doubling (sub-critical) and the 
exponent (3 with tangent-Hopf and the Double Hopf bifurcations. Such power-laws have 
been observed in various experimental studies of spatio-temporal intermittency |J. The 
values of the power-laws are listed in Table I. It is interesting to note that the values seen 
by us are in reasonable agreement with the power-laws seen in several experiments which 
involve quasi-one dimensional geometries fli~6 |, UTTf , ]Tj|. The exponent Ci which takes values 
between [1.9 — 2.2] for our system, is in good agreement with the laminar exponent seen in 



the case of Rayleigh-Benard convection in an annulus JT6|, whereas the exponent £3 which 
takes the values in the range [0.61 — 0.72] agrees well with the laminar exponent seen in 
the case of the roll coating system |T7|]. The exponent £ 2 which takes values in the range 
[1.3 — 1.35] and which is seen in the case of the periodic structure with long range spatial 
correlations is a lower bound on the laminar exponents seen in the case of the Rayleigh- 
Benard convection seen in a channel and in the Taylor Dean system where spatio-temporal 
intermittency is seen ptj, JT8]. Our numerical studies show that this is due to the fact that 



the temporally periodic structure with long-range spatial correlations undergoes a further 
bifurcation to spatio-temporal intermittency in the neighbourhood of e = 0.368,7 = 0-56 
for bifurcations from the synchronised fixed point x = 0.0. The distribution of laminar 
lengths in this region shows power-law behaviour with an exponent Cg ~ 1-63 which is in 
very good agreement with the laminar exponents observed for convection in a channel and 
for the Taylor-Dean system. This bifurcation does not appear to be local in nature. Details 
of this bifurcation are being explored further. 

All the results above are obtained for bifurcations from the synchronised fixed points 
for a lattice where the map parameter takes two distinct values one at each alternate site. 
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These results can be easily generalised to higher spatial and temporal periods and to lat- 
tices where the inhomogeneity has different periodicities. Bifurcations from higher spatial 
temporal periods can also give rise to spatio-temporal or spatial/temporal intermittency in 
many additional regions of parameter space. The question of whether the spatio-temporal 
intermittency with or without the accompanying power-laws for the distribution of laminar 
lengths persists in other regions of parameter space is being investigated further [|19|j . 

Thus we have shown that both spatial and spatio-temporal intermittency can arise in a 
inhomogeneous coupled map lattice. The phenomenon of intermittency is more widespread 
in inhomogeneous lattices than in the case of ordered lattices. The presence of pure spatial 
intermittency accompanied by temporally periodic behaviour is an interesting feature which 
arises in the case of the inhomogeneous system. The intermittency arises as a result of bifur- 
cations of co-dimension 2. Such bifurcations are also of interest in the case of other spatially 
extended systems ||15|| . The distributions of laminar lengths exhibit three distinct kinds of 
power laws, each associated with a distinct kind of bifurcation. The structure associated 
with the tangent-period-doubling bifurcation undergoes a further bifurcation associated with 
an additional power-law. The values obtained for these power-laws are in reasonable and 
intriguing agreement with those observed in a variety of experiments involving quasi-one 
dimensional geometries. We hence hope our analysis will be useful for the understanding of 
intermittent phenomena arising in other spatially extended systems, and in discussions of 
their genericity and universality. 
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TABLE 1.1 

Fixed Point x= 



Type 


Parameter region 


Eigenvalue 


P(l) « 


Nature 


Tangent Bifurcation 


_ (5 7 -21) (3 7 -9) 
(97-40) ' (77-24) 


+1 


Ci 


ST 


Tangent-PD Bifurcation 


e = 0.39,7 = 0.56 


+1, -1 


C2 


Spatial 


Tangent-Hopt Bifurcation 


e = 0.245,7 = 2.87 


a ± ib, +1 


c 3 


ST 



TABLE 1.2 

Fixed Point a: = 



Type 


Parameter region 


Eigenvalue 


P(l) w ^ 


Nature 


Tangent Bifurcation 


_ 57+8±(57- 1927+64) 2 
6 ~ 87 


+ 1 


Ci 


ST 


P D Bifurcation 


7 2 _3 7 _8±((7 2 -37-8) 2 -87(7-4) (7-6)) 3 


+ 1 


C2 


ST 


' ~~ (27(7-4)) 


Tangent - PD Bifurcation 


e = 0.63, 7= 1.180 


-1, +1 


C2 


Spatial 


Tangcnt-Hopt Bifurcation 


e = 0.06,7 = 2.9801 


+1, a±ib 


c 3 


ST 


Double-Hopf Bifurcation 


e = 0.5,7 = ^ 


a ± ib, c ± id 


c 3 


ST 



TABLE 2 



Experiment 


Value of laminar scaling exponent 


CML exponent £ 


Range of ( 


Convection in an annulus 


1.9 ±0.1 


Ci 


1.9-2.2 


Roll coating system 


0.63 ±0.02 


(3 


0.61-0.72 


Convection (channel) 


1.6 ±0.2 


Cg 


1.62 - 1.65 


Taylor Dean system 


1.67 ±0.14 


Cg 


1.62-1.65 
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I. FIGURE CAPTIONS 



1. Fig 1(a) shows the spatial intermittency in the lattice associated with tangent - 
period doubling bifurcation, i refers to position of site on the lattice. Fig 1(b) shows 
two iterates of the lattice. As can be seen both regular and irregular sites behave 
periodically in time. The fixed point is £ — 7e ~ M ~ 2 



2. Fig 2 The scaling behaviour of the laminar sites is shown. P(l) is the probability of 
obtaining a laminar region of length /. The length of the laminar region is defined as 
number of adjacent sites which remain within a particular accuracy. Three different 
exponents are found corresponding to three kinds of intermittency. The fixed point is 
f _ ye-fi-2 fj^g behaviour observed was obtained for a lattice of size 50,000 iterated 

^ fi— ■ye 

for 20000 iterates starting with 100 random initial conditions. The accuracy to which 
the laminarity of the region was checked was 10 -5 . The asterisks denote the behaviour 
with exponent Ci, crosses C2 and pluses denote £3- The ranges of the exponents are 
given in the caption of Table I. £3 shows departures from this power-law over the third 
decade. The axes are marked in the natural log-scale. 

3. Fig 3 The spatial second return map x(i + 2) vs x(i) where i is the site index is plotted 



for £ 



II. TABLE CAPTIONS 



1. Table 1.1 We list the bifurcations from the synchronised fixed point x = 0. The type 
of bifurcation involved, the region in parameter space where the bifurcation conditions 
are satisfied, the manner in which the eigenvalue crosses unit circle and the power law 
exponent P(l) fa fa are listed. The range of the exponent £1 is [1.9 — 2.2], (2 lies in 
the range [1.3 — 1.35] and (3 in the range [0.61 — 0.72]. We also identify the nature 
of the intermittency, whether spatial or spatio-temporal(ST). In Table 1.2 the same 
quantities are given for the other fixed point. 

2. Table 2 This table lists the values of the spatial laminar exponent observed in ex- 
periments (values as quoted in Ref. f|) and compares their values with the laminar 
exponents of our CML model. It is clear that the exponents £1 and £3 are directly ob- 
served in fluid experiments on convection in an annulus and in the roll coating system. 
As explained in the text, the exponent (2 serves as a lower bound on the ex ponent 
observed in convection in a channel and in the Taylor Dean system. The experimen- 
tally observed exponent actually coincides with the CML exponent (g which occurs 
when the structure associated with the exponent (2 undergoes a further bifurcation. 
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Fig. 1(a) AS and NG 
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